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IONIZATION  AND  PULSE  LETHARGY  EFFECTS 
IN  INVERSE  CHERENKOV  ACCELERATORS 

1.  Introduction 

Although  high  power  sources  of  electromagnetic  radiation,  such  as  lasers,  are  capable  of 
proving  extraordinarily  high  electric  fields  for  the  acceleradon  of  particles,  numerous 
fundamental  and  technological  issues  must  be  resolved  before  practical  high  gradient  accelerators 
can  be  realized  [1-1 1].  There  are  three  fundamental  issues  tiiat  must  be  addressed  in  any  type  of 
high  gradient  accelerator.  These  are:  i)  radiation  beam  guiding  over  extended  distances,  ii)  phase 
coherence  (phase  velocity  slippage)  between  die  accelerated  particles  and  the  electromagnetic 
field,  and  iii)  material  ionization  arising  from  either  the  driver  fields  or  from  the  accelerated 
particles. 

In  this  paper  we  analyze  and  discuss  the  fundamental  issues  as  tiiey  relate  to  the  inverse 
Cherenkov  accelerator  (ICA)  [12-18].  The  electromagnetic  driver  in  the  ICA  can  be  an  intense 
laser  or  millimeter  wave  pulse.  In  free  space  and  in  the  absence  of  boundaries  the  phase  velocity 
of  the  accelerating  field  exceeds  the  vacuum  speed  of  light,  c.  As  a  consequence,  particles 
continually  phase  slip  relative  to  the  field,  Le.,  phase  velocity  slippage  and,  after  a  slippage 
distance,  acceleration  ceases.  In  the  ICA  the  phase  velocity  of  the  accelerating  field  is  reduced  by 
introducing  a  gas  in  the  accelerating  region  or  by  lining  the  interior  of  a  waveguide  with  a 
dielectric  material 

A  dielectrically  lined  or  gas-filled  waveguide  ICA  configuration,  i)  avoids  diffraction  of  the 
driving  electromagnetic  beam,  ii)  overcomes  electron  slippage,  and  iii)  the  acceleration  in  the 
dielectrically  lined  ICA  is  in  a  vacuum.  However,  the  large  electric  fields  associated  with  short, 
intense  electromagnetic  pulses  can  readily  ionize  the  dielectric  material  or  gas  [5,19-22].  Partial 
volume  ionization  will  modify  the  dispersive  characteristics  of  the  waveguide,  due  to  plasma 
Manuscript  approved  December  12,  1996 
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formation.  A  small  amount  of  ionization  can  result  in  the  disruption  of  the  acceleration  process 
since  the  phase  velocity  of  the  clectroma^etic  fields  will  be  altered,  resulting  in  phase  velocity 
slippage.  Ionization  effects  therefore  impose  an  upper  limit  on  the  pulse  duration.  However,  the 
driving  pul^  cannot  be  arbitrarily  short  since  the  pulse,  having  a  group  velocity  less  than  c,  can 
slip  behind  the  accelerated  electrons.  This  effect  is  referred  to  as  group  velocity  slippage  or  pulse 
lethargy  and  places  a  loww  linut  on  the  pulse  duration.  The  combination  of  ionization  and  pulse 
lethargy  effects  impose  severe  limitations  on  the  energy  gain  in  inverse  Cherenkov  accelerators. 
For  sufficiently  high  current  densities,  the  self-fields  of  the  accelerated  electron  beam  may  also 
result  in  ionization.  Self-field  ionization  processes  will  not  be  considered  in  this  paper. 

This  paper  is  organized  as  follows.  In  Section  n  we  dwive  and  discuss  various  limits  of 
the  dispersion  relation  for  transverse-magnetic  (TM)  modes  of  a  dielectrically-lined  optical 
waveguide.  This  TM  mode  consists  of  a  radial  and  an  axial  electric  field  as  well  as  an  a^.uthal 
magnetic  field.  The  axial  electric  field,  which  is  responsible  for  acceleration,  is  peaked  on  axis, 
while  the  radial  electric  field  vanishes  on  axis.  In  Section  ffl  we  derive  the  critical  average  plasma 
density.  When  the  average  plasma  density  reaches  this  critical  value  the  dispersive  properties  are 
modified  and  the  phase  velocity  of  the  electromagnetic  wave  is  significantly  altered,  resulting  in 
phase  velocity  slippage.  In  Section  IV  the  buildup  of  plasma  density  within  the  dielectric  material 
is  analyzed.  Here  we  assume  that  ionization  takes  place  within  the  dielectric,  Le.,  surface  effects 
are  neglected.  This  analysis  is  based  on  a  rate  equation  for  plasma  generation,  which  includes 
tunneling  and  collisional  (avalanche)  ionization.  In  Section  V  we  obtain  an  expression  for  the 
limit  on  the  electromagnetic  pulse  duration  in  terms  of  the  accelerating  gradient  The  case  of  a 
gas-filled  waveguide  ICA  is  analyzed  in  Section  VI.  The  maximum  accelerated  elecuon  energy 
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imposed  by  group  velocity  slippage,  Le.,  pulse  lethargy  and  ionization,  is  obtained  in  Section  VIL 
Section’ Vin  presents  numerical  results  for  the  pulse  duration  and  accelerating  gradient,  limited  by 
ionization  and  pulse  lethargy,  for  a  10  |Jtm,  1  mm  and  1  cm  wavelength  driver  in  the  two  ICA 
configurations.  The  summary  and  discussion  is  given  in  Section  DC. 

n.  Fiftlds  and  Dispersion  Relation  for  Inverse  Cherenkov  Accelerator  vnlh  Dielectric  Liner 
A.  Axially  Symmetric  Transverse  Magnetic  Mode 

The  dispersion  characteristics  of  a  dielectrically  lined  waveguide  is  analyzed  in  this  section. 
A  dielectrically  lined  waveguide  can  support  an  axially  symmetric  transverse-magnetic  (TM) 
mode.  The  cross-sectional  view  of  the  dielectrically-lined  waveguide  is  shown  in  Rg.  1,  where 
the  inner  surface  of  the  dielectric  is  at  r  =  a  and  the  conducting  outer  surface  is  at  r  =  b.  The 
particular  TM  mode  under  consideration  consists  of  radial  and  axial  electric  fields  as  well  as  an 
azimuthal  magnetic  field.  The  field  components  within  the  central  vacuum  region  (0  <  r  <  a)  are 


E,  =  Eo(Ji(kxr)/ Ji(kxa)X(z,t)+c.c.,  (la) 

=i(kj./k)Eo(Jo(kxr)/ Ji(kxa)X(z.t)+c.c.,  (lb) 

Be  =(o)/ck)Eo(Ji(kj.r)/ Ji(kj^a)X(z,t)-i-c.c.,  (Ic) 

and  within  the  outer  dielectric  region  (a  ^  r  <  b), 

Ej  =  (AJ,((xr)+BY,(otr)X(z.t)  +  c.c.,  (2a) 

Ej  =  i(a  /  k)(AJo(ar) + BYo(ar)X(z.  t) + c.c.,  (2b) 

Be  =e(a>/ck)(AJ|(ar)+BYj(ar)X(z»t)-»-c.c.,  (2c) 


where  ©  is  the  frequency,  k  is  the  axial  wavenumber,  f(z,t)  =  (l/2)exp(i(kz-  ©t)), 
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kj^  =  (CD^  /  is  the  transverse  wavenumber  within  the  inner  region,  a  =  (eco^/c^  —  k^)‘^ 

is  the  transverse  wavenumber  within  the  outer  region,  e  is  the  dielectric  constant.  Eo,  A,  B  are  the 
constant  amplitudes,  Jb(x)  and  Yb(x)  are  the  ordinary  Bessel  functions  of  order  n  and  c.c.  denotes 
the  complex  conjugate.  The  accelerating  axial  electric  field  given  by  Eq.  (lb)  is  peaked  along  the 
z-axis. 

For  application  to  electron  acceleration  we  will  consider  the  regime  where  the  phase 
velocity  is  near  c,  i.e.,  Yph  s  ea  /  k  =  c.  For  phase  velocities  equal  to  c,  i.e.,  kj. = 0,  the  axial 

electric  field  is  independent  of  radial  position  and  given  by 


E,=i— Eof(z,t)+c.c.,  (3) 

Tta 

where  X  =  27tc/o)  is  the  wavelength.  The  axial  and  radial  components  of  the  electric  field  for 
phase  velocity  equal  to  c  are  shown  in  Fig.  2. 


B.  Dispersion  Relation 

Applying  the  appropriate  boundary  conditions  at  r  =  a  and  r  =  b  we  obtain  the  following 
dispersion  relation  for  the  TM  mode. 


JoC^^a 
^  J,(kj.a)  4 


where  a^=eci)^/c^-k^=kj^+(e-l)o)^/c*  and 


(4) 


4a  Y  Jo  (cxb)Yo(oca)  -  Jo(aa)Yo(ocb)^ 
ea>\  Jo(ab)Y,((xa)-Ji(aa)Yo(ab) ; 


In  addition,  the  amplitudes  A  and  B  in  terms  of  Eo  are 
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(5a) 


_ Yo(ab) _ 

e  Jo(ab)Yi(aa)-Ji(aa)Yo(ab)’ 

B  =  -AJo(cxb)/Yo(ab).  (5b) 

For  phase  velocities  close  to  c,  the  dispersion  relation  can  be  simplified  by  taking  k^a  to  be  small 

on  the  left-hand  side  of  Eq.  (4),  giving 

0)2  /  c^  -  k^  =  8 /  -  F(ao).  (6) 


where  tto  =  (e  - 1)'^^0)  /  c .  In  obtaining  Eq.  (6)  the  small  argument  forms  Jo(x)  =  1  -  x^M  and 
Ji(x)  =  x/2  -  xVl6  were  used. 

The  dispersion  relation  in  Eq.  (6),  which  is  valid  for  phase  velocities  near  c,  can  be  further 
simplified  by  considering  several  limiting  regimes.  A  relevant  regime  corresponds  to  the  limit 
where  the  waveguide  is  large  and/or  the  dielectric  constant  is  far  from  unity,  corresponding  to  the 
inequality  ttQa  »  1,  i.e., 

X/a«27i(e-l)*'2.  (7) 

The  dispersion  relation  in  this  regime  is  given  by 

fl)2/c2-k2  =8/a2-^^tan(ao(b-a)).  (8) 

ea  '  ' 

In  obtaining  Eq.  (8)  the  asymptotic  form  F(ao)  =  (4ao  /  ea)  tan(ao(b  —  a))  has  been  used.  The 


components  of  the  electric  field  within  the  dielectric  are 

cos[ao(b-r)] 


e  • - - ^f(z,t)  +  C.C., 

cos[ao(b-a)] 


iX  an[ao(b-r)] 


[na  an[ao(b  -  a)) 


f(z,t)+c.c. 


(9a) 

(9b) 
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The  dispersion  relation  in  Eq.  (8)  can  be  solved  explicitly  in  three  cases.  To  obtain  the 
frequencies  associated  with  these  cases  it  is  convenient  to  rewrite  the  dispersion  relation  in  Eq. 
(8),  for  phase  velocities  equal  to  c,  in  the  form 

.  tan|  =  2e(b-a)/(a^).  (10) 

where  ^  =  ao(b  -  a).  In  Case  1  the  right-hand  side  (RHS)  of  Eq.  (10)  is  much  greater  than  unity, 
i.e.,  X/a  » 7t(e  -  l)*^/e,  and  ^  «=  (/+ 1  /  2)7C,  where  t  =  0,1,2 .  This  case  is  compatible  with 


the  inequality  in  Eq.  (7)  provided  e  is  sufficiently  large  compared  to  unity.  In  Case  2  the  RHS  is 
small  compared  to  unity,  Le.,  X/a  « Jc(£  - 1)*^/£,  and  ^  *  frt,  where  /  is  a  large  integer.  Note 
that  the  inequality  defining  this  case  is  more  restrictive  than  that  in  Eq.  (7).  Finally,  Case  3 
corresponds  to  %  «  1,  i.e„  when  the  liner  is  veiy  thin,  (b  -  a)/X  « (e  -  1)  '^/(27t).  We  emphasize 
that  the  inequality  in  Eq.  (7)  applies  to  all  three  cases.  The  wavelengths  associated  with  these 
cases  are  obtained  from  Eq.  (10), 


(lla) 

(llb) 

(llc) 


where  Eq.  (1  la)  corresponds  to  Case  1  with  t  nodes  in  the  liner,  Eq.  (1  lb)  corresponds  to  Case 
2  with  t  nodes  in  the  liner,  and  Eq.  (1  Ic)  corresponds  to  Case  3  with  no  nodes  in  the  liner.  The 


radial  component  of  the  electric  field  inside  the  dielectric  material  is  given  by 

(-l)^(X  /  jca)(e  /  (e  -  +  c.  c. , 


(-1)^  cos^l  ^^jf(2,  t)  +  C.  C. , 
f(z,t)  +  c.c.. 


(12a) 

(12b) 

(12c) 
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where  ^  in  Eqs.  (12a)  and  (I2b)  is  given  by  Eq.  (10)  in  the  appropriate  limit  Similarly,  the  axial 
components  of  the  electric  field  inside  the  dielectric  liner  is  given  by 


Ej  *=E, 


ArJ  jca 


(-1)^  an^^^^jf(z,  t)  +  C.C., 

(-l)^(jca  /  X)((e  - 1)’^^  /  e)an^^-^^^(z,  t) + c.c.. 


b-r 


b-a 


f(2,t)  +  C.C. 


(13a) 

(13b) 

(13c) 


Using  these  expressions  one  can  readily  verify  that  the  field  components  satisfy  the  appropriate 
continuity  and  boundary  conditions  at  the  surface  of  the  liner  (r  =  a)  and  at  the  conducting  wall  (r 


=  b). 

Other  regimes  of  the  dispersion  relation  in  Eq.  (6)  can  be  evaluated  but  do  not  lead  to 
phase  velocities  equal  to  c.  Simplification  of  the  dispersion  relation  in  Eq.  (6)  in  the  regime  of 
long  wavelengths  may  be  made  by  taking  Oob  «  I,  i.e.,  X/b  »  2jt(£  - 1)’^.  This  regime  is 
incompatible  with  phase  velocities  equal  to  c,  as  is  shown  below.  The  dispersion  relation,  in  this 
regime,  takes  the  form 

/  c^  -  k^  =  8  /  a^  -  E"‘  (2ao  f  logc(b  /  a). 


The  solution  of  this  dispersion  relation  >n1ien  the  phase  velocity  is  equal  to  c  can  be  written  in  the 
form 

To  be  consistent  with  the  inequality  defining  this  regime  we  require  X  /  (27ib(E  - 1)^^^  )  » 1, 

which  is  impossible  since  b  >  a  and  e  >  1.  Hence,  there  is  no  solution  to  the  dispersion  relation  as 
X  becomes  large  compared  to  the  dimensions  of  the  outer  conductor.  The  final  regime  tiiat  we 
identify  explicitly  corresponds  to  Ooa  «  1  and  Oob  »  1.  Again  this  regime  is  incompatible  with 
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phase  velocities  equal  to  c  since  the  dispersion  relation  reduces  to  2£  —  -  (otoa)  iogeCotoa).  In  the 
limit  oto  0,  the  RHS  vanishes  and  hence  no  solution  to  the  dispersion  relation  exists. 


in.  Critical  Average  Plasma  Density 

The  large  fields  needed  for  accel^^don  can  ionize  the  dielectric  material  and  change  the 
dispersive  properties  of  the  waveguide.  If  die  dispersive  properties  change  sufficiently  phase 
velocity  slippage  will  occur  and  the  acceleration  will  cease.  As  ionization  takes  place  the 

dielectric  constant  changes  from  e  to  e  +  8e,  where  5e  =  -  to J  /  (D^  is  the  contribution  to  the 

i/2 

dielectric  constant  due  to  ionization,  tOp  =  ^4jtq^^np^  /  m  j  is  the  average  plasma  frequency,  q 
is  the  charge  and  m  is  the  mass  of  an  electron  and  (np)  is  the  average  plasma  density  within  the 
waveguide,  i.e.,  ts  the  peak  plasma  density  times  the  ratio  of  the  cross-sectional  area  of  the 

plasma  to  the  cross-sectional  area  of  the  waveguide.  When  the  change  in  F(o(o)  in  Eq.  (6)  due  to 
ionization  equals  ~  8/a^  the  optical  properties  of  the  guide  change  significantly  and  phase  velocity 
slippage  disrupts  the  acceleration.  Ionization  in  the  dielectric  material  will  not  disrupt  the 
acceleration  provided 

|5e|  =  ©5  /  CD^  « (8  /  a^)|aF  /  ae|”*.  (14) 

This  inequality  defines  a  critical  average  plasma  density  within  the  dielectric,  iierit,  given  by 

where  r*  =  qVmc*  =  2.8  x  10*’’  cm  is  the  classical  electron  radius.  If  the  average  plasma  density 
^np^  generated  within  the  dielectric  exceeds  ticnt  the  acceleration  is  disrupted  because  of  phase 

velocity  slippage. 
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The  critical  average  plasma  density  can  be  evaluated  in  the  limiting  cases  discussed  in 
connection  with  the  dispersion  relation  in  Section  IL  The  critical  average  plasma  density  in  the 
duee  cases  is 


(16a) 

(16b) 

(16c) 


where  Eqs.  (16a)-(16c)  refer  to  C!ase  1-3,  respectively.  As  an  illustration  we  take  e  =  2,  a  =  1(X) 
urn,  a/(b-a)  =  2  and  find  that  the  critical  density  for  phase  slippage  in  (Zase  2  is  ncnt  =  5  x  10*‘ 
on'^  and  is  independent  of  the  wavelength.  For  a  10  |im  wavelength,  this  value  is  more  than  two 
orders  of  magnitude  smaller  than  the  plasma  density  for  complete  reflection  (~10”  cm’’). 


IV.  Ionization  in  Solids 

The  rate  of  change  of  plasma  density  (density  of  free  electrons)  is  given  by 

^  =  S  +  Wnp-v,np,  (17) 

where  S  is  the  photo-ionization  source  due  to  electron  tunneling,  W  is  the  collisional  ionization 
rate,  and  v,  is  the  recombination  rate.  The  photo-ionization  source  acts  as  an  initial  source  of  free 
electrons  which  are  further  increased  in  number  by  collisional  processes,  i.e.,  by  electron 
avalanche.  Typically,  for  parameters  of  interest  here,  the  recombination  rate  is  small  compared  to 
the  collisional  ionizadon  rate,  i.e.,  Vr  <  W,  and  will  be  neglected.  Note  that  S  and  W  have  units  of 
sec*’  cm*’  and  sec*’  respectively  and  will  be  specified  later.  The  solution  of  Eq.  (17)  is 
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Hp  =  (S  /  W)(exp(Wt)  - 1)+  Hpo  exp(Wt), 


(18) 


where  Dpo  is  the  seed  electron  density  and  it  is  assumed  that  die  plasma  density  remains  small 
compared  to  the  neutral  density. 

A.  Collisional  Ionization 

The  collisional  ionization  rate  can  be  estimated  by  using  a  classical  free  electron  model  for 
the  electron  energy  gain.  The  rate  of  change  of  energy  of  an  electron  undergoing  collisions  in  the 
presence  of  a  temporally  periodic  electric  field  is  given  by  [21^] 

(19) 

dt  2m(0)^  +  Vm)  M 

where  U  is  the  electron  energy,  E  is  the  peak  electric  field  amplitude,  v„  is  the  momentum 
transfer  frequency,  and  M  is  the  mass  of  the  neutral  atoms.  The  first  term  on  the  right-hand  side 
of  Eq.  (19)  represents  the  rate  of  energy  increase  of  an  electron  undergoing  collisions  while  being 
accelerated  in  the  alternating  field.  The  second  term  represents  the  energy  damping  due  to 
inelastic  collisions  between  electrons  and  neutral  atoms.  Note  that  E  in  Eq.  (19)  refers  to  the 
peak  electric  field  in  the  dielectric.  Equation  (19)  indicates  that  the  electron  energy  increases 
monotonically  and  saturates  because  of  inelastic  electron  collisions  with  neutral  atoms.  Solving 
Eq.  (19),  we  obtain 

U(t)  =  Usat  (1  -  exp(-Vot)) ,  (20) 

where  U,at  =  (v^  /  Vo)(l + vj,  /  a)^)~*Uos  is  the  saturation  (maximum)  electron  energy, 

Uqs  =  nivos  /  2  is  the  electron  oscillation  energy,  Vos  =  qE/m(0  is  the  oscillation  velocity  and  Vo  = 


10 


2mVn/M  is  the  energy  damping  frequency.  Typically,  in  solids  Vn  ~  10‘*  sec’*  [23]  and  Vo  ~  10“ 
sec'*.  When  the  electron  energy  reaches  the  ionization  wiergy  Ui,  die  dectron  collisionally  ionizes 
the  atom.  The  time  for  this  to  occur  is  the  ionization  time  W*.  Setting  U  =  Ui  and  t  =  W*  in  Eq. 
(20)  determines  W, 

W  =  -Vo/loge(l-U,/U^),  (21) 

for  UiAJ«  <  1  and  W  =  0  for  Ui/U,^  >  1.  For  Ui/U«< «  1.  Eq.  (21)  becomes 
W  =  (Uos  /  Ui  )Vn,  /  (vj,  /  0)2  + 1).  Note  that  (U,  AJo.)*'^  =  it  is  the  Keldysh  parameter  [24]  and 

•yi  <  1  corresponds  to  die  tunneling  ionization  regime.  As  noted  earlier,  we  assume  that  the 
ionization  takes  place  within  the  dielectric  material  The  minimum  electric  field  for  ionization, 
associated  with  long  pulses,  can  be  estimated  from  Eq.  (20)  by  taking  U(t  -¥  <•«•)  =  Ui: 

E  =  (2m  /  M)’'2[(co2  +  (22)  . 

where  Ui  =  Ui/Uh  is  the  atomic  ionization  enei;gy  Ui  normalized  to  that  of  hydrogen, 

Uh  =  /  (2ao )  =  13.6  eV,  E  =  E  /  Eji  is  the  peak  electric  field  normalized  to  the  hydrogenic 

electric  field  En  =  |q|  /  ao  =  5.2  GV  /  cm,  Gq  =  ^jcUh  /  h  =  4.1  x  10*^  sec”*  is  the  atomic 
frequency,  Uq  =  (h  /  27tq)^  /  m  =  5.3  x  10“®  cm  is  the  Bohr  radius,  and  h  is  Planck’s  constant 

Irregularities  on  the  surface  of  the  dielectric  liner  may  permit  surface  plasma  formation  at  electric 
field  amplitudes  substantially  smaller  than  dial  implied  by  Eq.  (22).  The  results,  in  this  paper,  for 
maximum  accelerating  fields  are  therefore  expected  to  be  upper  bound  estimates.  The  rough 
estimate  for  the  coUisional  ionization  rate  in  Eq.  (21)  is  used,  together  with  the  tunneling  source 
terra  S,  to  obtain  the  plasma  density  in  Eq.  (18). 
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B.  Tunneling  Ionization  Source 

The  ionization  source  in  solids,  averaged  over  times  much  longer  than  the  field  period 
27t/ciX  in  the  tunneling  regime  7k  =  (UiAJoi)*^  <  1,  is  given  by  [25] 


SsA,E^^^exp[-p,/E], 

(23) 

where 

A,=L5xl0^/U,*'^ 

(24a) 

p,  =0.56Ui^'2. 

(24b) 

In  Eqs.  (23)  and  (24),  A*  has  units  of  sec**  cm*’' 

Figure  3  shows  the  tunneling  ionization  source  tern  and  the  collisional  ionization  rate  as 
functions  of  the  normalized  electric  field.  For  pulses  that  are  shorter  than  a  few  picoseconds, 

collisional  ionization  is  negligible  for  W  <  10***  sec**.  The  precipitous  drop  in  W  below  E  =  2x10 
occurs  as  Umax  approaches  Ui  in  Eq.  (21).  The  tunneling  terra  S  does  not  lead  to  significant 

ionization  for  E  <  5xl0'^  but  ionization  increases  rapidly  for  higher  fields. 


V.  Limitation  on  Laser  Pulse  Duration  and  Accelerating  Gradient 

To  avoid  modifying  the  dispersive  properties  of  the  waveguide  the  average  induced 
plasma  density  must  be  less  than  the  critical  average  density  ncrit»  i-©-  (*^p}  « ‘'crit-  Using 
Eq.  (18)  for  np  we  find  that 
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(np)  =  ((S  /  W)(exp(WtL)  - 1)}+ npo(exp(WTL))  «  nent. 


(25) 


where  (  )  denotes  the  average  over  the  cross-sectional  area  of  the  waveguide  and  Tl  is  the  laser 

pulse  duration.  The  average  plasma  density  in  Eq.  (25)  consists  of  a  tunneling/collisional 
ionization  term  and  a  contribution  from  the  seed  electron  density.  To  perform  the  cross-sectional 
average  we  note  that  in  the  ionization  part,  the  source  function  S  is  a  highly  sensitive  function  of 
the  field  amplitude  through  the  exponential  dependence  on  the  field  [see  Eq,  (23)].  Because  of 
this  sensitivity,  the  spatial  average  of  the  ionization  term  in  Eq.  (25)  is  applied  only  to  the  function 
S  while  the  collisional  ionization  rate  W  is  evaluated  at  the  peak  field  amplitude.  The  cross- 
sectional  average  of  the  contribution  from  the  initial  plasma  electron  density  can  be  written  as 
(exp(WxL))  =  f  exp(WjjjtL)  where  f  is  a  filling  factor  due  to  the  spatial  variation  of  the  field 
within  the  dielectric.  Employing  these  rq)proximations,  Eq.  (25)  becomes 

w;'(«p(W„Tl)-1XS>+  fhpo  exp(W„TL) «  n^,,  (26) 


where  Wm  is  the  maximum  value  of  W,  Le.,  W  evaluated  at  die  peak  electric  field  Em  inside  the 
dielectric.  The  filling  factor  f  is  found  to  be  given  by 


f  = 


1/2, 

2(b-a)/b, 


(27a) 

(27b) 

(27c) 


where  Eqs.  (27a)  refers  to  either  Case  1  or  Case  2  when  WbXl  «  1,  Eq,  (27b)  refers  to  either 
Case  1  or  Case  2  when  WbXl  »  1,  and  Eq.  (27c)  refers  to  only  Case  3. 
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The  pulse  duration  of  the  driver  in  the  ICA  roust  not  exceed  a  critical  pulse  duration  Xcni- 
Substituting  Eq.  (16)  into  Eq.  (26)  we  obtain  the  following  inequality  on  the  pulse  duration 

«  ’^erit  =  logc[('CoWn,  +  !)/(!  +  fopoWn,  /  (S))]  ,  (28) 

where  to  =  nethf<S>.  If  the  laser  pulse  duration  is  small  compared  to  the  collisional  ionization 
time,  Le.,  Tl  «  W“* ,  and  the  initial  plasma  density  is  zero,  Upo  =  0,  the  inequality  in  Eq.  (28) 
becomes  Xl  «  Tq  =  n^nt  /  (S) .  The  inequality  in  Eq.  (28)  places  a  limit  on  the  pulse  duration, 
which  is  a  function  of  the  accelerating  gradient,  E^  =  (X  /  a7c)Eo .  If  the  pulse  duration  exceeds 

the  value  given  by  Eq.  (28),  sufficient  ionization  will  take  place  to  modify  the  dispersive 
properties  of  the  waveguide,  causing  phase  velocity  slippage  to  terminate  the  acceleration. 

The  critical  time  tent  in  Eq.  (28)  must  be  determined  explicitly  for  the  three  cases.  The 
average  ionization  source  terra  is 

b 

(S)  =  27iJ  Srdr  /  (nh^ ),  (29) 

a 

where  the  radial  variation  in  S  arises  from  that  of  the  electric  field.  Since  the  ionization  source 
terra  S  is  a  highly  sensitive  function  of  the  electric  field  due  to  the  exponential  factor  in  Eq.  (23), 
the  integral  in  Eq.  (29)  can  be  approximately  evaluated.  In  the  limit  defined  by  the  inequality  in 
Eq.  (7)  the  electric  field  is  given  by  Eqs.  (12)  and  (13).  Making  use  of  these  expressions  in  Eq. 
(29)  and  expanding  about  each  of  the  £  peaks  for  Cases  1  or  2,  and  about  the  peak  at  r  =  a  for 
Case  3,  we  obtain 
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(30a) 


<S)  = 


jr 

2b 


2 


8E„/p,. 


(30b) 


where  Eq.  (30a)  refers  to  either  Case  I  or  Case  2  and  Eq.  (30b)  refers  to  only  Case  3  and  S„  is  S 


evaluated  at  E„,  =  E^,  /  E^,  Le., 
Sm=A,En,^'^expj 


5/2 - ^[-pj/Eni]. 


(31) 


Using  Eq.  (30)  for  (S)  and  Eq.  (18)  for  Ucnt.  we  obtain 

v3 


'Co  =  ncrit^(S)  = 


e  1  a 


8n  a%  b-aa'^Sj 


(32a) 

(32b) 


b-a' 

lPs/E„.  (32c) 

where  Eq.  (32a)  refers  to  Case  1,  Eq,  (32b)  refers  to  Case  2,  and  Eq.  (32c)  refers  to  Czse  3. 
From  Eq.  (28),  values  of  Xo  for  the  various  cases  can  be  used  to  determine  the  limits  on  the  pulse 
duration  and  the  accelerating  gradient 


VI.  Gas-Filled  Inverse  Cherenkov  Accelerator 

Another  configuration  for  the  inverse  Cherenkov  accelerator  consists  of  a  waveguide  that 
is  filled  with  a  neutral  gas.  The  presence  of  the  neutral  gas  modifies  the  electromagnetic 
properties  of  the  waveguide  making  the  pha^  velocity  of  the  accelerating  electric  field  equal  to 
the  speed  of  light  For  sufficiently  low  gas  densities,  the  collisional  scattering  of  the  accelerated 
electrons  can  be  neglected  [26]. 
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In  this  section  we  consider  the  limitations  on  electromagnetic  pulse  duration  and 
accelerating  gradient  of  an  ICA  consisting  of  a  gas-filled  cylindrical  waveguide  of  radius  a.  For 

the  TM  mode  the  fields  are  given  by  Eq.  (1),  with  k  /  c^  -  ,  where  e  >  1  is  the 

dielectric  constant  of  the  gas.  The  dispersion  relation  for  the  TM  mode  is  given  by  kj^a  =  pon.  or 

£©^  / c*  -k^  =  (Pon  / a)^» 

where  pod  is  the  nth  zero  of  Jo.  Writing  e  =  1 + Ae,  where  Ae  >  0  is  the  contribution  of  the  neutral 
gas  to  the  dielectric  constant,  from  Eq.  (33)  we  find  that,  for  the  phase  velocity  to  equal  c,  the 
wavelength  is  given  by 

X/a  =  (27t/pon)Ae‘'^  (34) 

Since  Ae  «  1  for  gases  at  moderate  pressures,  i.e.,  a  few  atmosphere,  the  expression  in  Eq.  (34) 
implies  that  X/a  «  1.  Noting  that  E*  =  (X/Jca)Eo  on  axis,  the  accelerating  field  is  necessarily  small 
compared  to  the  field  amplitude  Eo. 

To  determine  the  critical  average  plasma  density  ncnt  in  the  gas-filled  ICA,  we  proceed  as 

2  2 

in  Sec.  EH.  Accounting  for  the  ionization  of  the  ambient  gas,  we  write  e=l+AE-  ©p/©  and 
obtain  the  critical  plasma  density 

=  jcAe  /  (r^X? )  =  Pon  /  47ta^rg.  (35) 

Typically  Ae  ~  lO"*  for  a  gas  at  atmospheric  pressure,  and  for  a  X  =  10  iim  laser-driven  gas-filled 
ICA,  the  critical  density  is  found  to  be  ncnt  =  lO’^  cm’. 

The  development  of  plasma  in  the  gas,  undergoing  ionization  due  to  tunneling  and 
coUisional  processes,  is  given  by  Eq.  (17).  The  ionization  source  term  in  Eq.  (17)  for  a  gas  is 
given  by  [24] 
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(36) 


S  =  AgE"^'^exp(-pg/E). 

where 

Ag=L6xl0”u/'V„0. 

pg=0.67Up.  (37b) 

Dgo  is  the  neutral  gas  density  in  units  of  cni*^  and  \  in  Eq.  (37a)  is  in  units  of  sec  ’  cm  .  The 
collisional  ionization  rate  W  is  given  by  Eq.  (21).  Similarly,  the  rate  of  change  of  electron  energy 
is  given  by  an  Eq.  (19),  where,  the  momentum  transfer  frequency  v„  in  a  gas  is  proportional  to  the 
pressure  and  is  on  the  order  of  10*^  sec*'  at  atmospheric  pressure  [21,22].  The  limitation  of  the 
pulse  duration  and  accelerating  gradient  in  the  gas-filled  ICA  is  given  by  Eq.  (28).  For  a  gas-filled 
ICA,  to  =  ncrit/<S>,  where  ricnt  is  given  by  Eq.  (35),  <S>  is 

<S)  =  S„(2JlE„/Pg)‘'^  (38) 


and  Sm  is  the  value  of  S  evaluated  at  peak  electric  field  Eo,  i.e.,  Eq  =  E^,  /  Ey.  The  filling  factor 
f  in  Eq.  (28)  for  the  gas-filled  waveguide  ICA  is 


f 


1/2, 

(Jt/2)''^(W„tL)-''^ 


(39) 


where  the  upper  (lower)  value  is  for  W^tL  « 1  (» 1). 


Vn.  Group  Velocity  SUppaee  (Pulse  Lethargy) 

Ionization  effects  place  an  upper  limit  on  the  pulse  duration  of  the  electromagnetic  driver, 
see  Eq.  (28).  Group  velocity  slippage  or  pulse  lethargy,  on  the  other  hand,  places  a  lower  limit  on 
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the  pulse  duration  of  the  driver.  That  is,  the  distance  the  pulse  slips  behind  the  accelerated 
electrons  must  be  less  tiian  the  pulse  length, 

where  Vj  is  the  group  velocity  and  L  is  die  propagation  distance.  Using  Eq.  (8)  the  group  velocity 
Vg  s  3©  /  3k  is  given  by 

l-f  (X  /  4;i)(e  -  1)~*^^3F  /  3ao 
1 + (X  /  4jt)e(e  -  ir‘'^3F  /  3ao  ’ 

where  F  =  (4a</ea)tan(ao(b  -  a)).  Making  use  of  Eq.  (41),  the  group  velocity  in  the  three  cases  is 


given  by 

1/e. 

l  +  2(b-a)/(ea) 
''.''“I  l+2(b-a)/a  ■ 
l  +  4(b-a)/(ea) 

■  l  +  4(b-a)/a 


(42a) 

(42b) 

(42c) 


Note  that  Eq.  (42a) ,  which  corresponds  to  Case  1,  indicates  that  the  group  velocity  »  c/e.  This  is 
consistent  with  the  relationship  VphV*  =  cVe  for  a  completely  fiUed  waveguide.  Substituting  Eq. 
(42)  into  Eq.  (40)  yields  the  limit  on  the  interaction  length  due  to  pulse  lethargy, 

L<G(e,a,b)ctL,  (^3) 


where 

1, 

l+ea/2(b-a), 
l+ea/4(b-a). 


G  =  (e-1)‘^ 


(44a) 

(44b) 

(44c) 
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In  obtaining  Eq.  (41)  we  have  neglected  the  effects  of  the  induced  plasraa,  Le.,  tq,  —  0.  For  a  gas- 
fiUed  ICA,  discussed  in  Section  VI,  the  group  velocity,  neglecting  transverse  effects,  is  given  by 
Vg  /  c  £~^,  which,  upon  substitution  into  Eq.  (40)  yields  the  linut  on  the  interaction  length, 

L<(e-l)-’cXL,  (45) 

which  is  identical  in  form  to  Case  1,  Eq.  (44a),  of  the  dielectrically-lined  configuration. 

The  limitation  on  the  pulse  length  due  to  ionization,  Eq.  (28),  together  with  the  limitations 
imposed  by  pulse  lethargy,  Eqs.  (43),  present  a  severe  constraint  on  the  interaction  length,  and 
therefore  the  maximum  energy  gain,  in  die  ICA  The  limitadon  on  die  interaction  length  L  can  be 
stated  in  terras  of  the  maximum  electron  energy  Ubbx  and  the  acceleration  gradient  E^.  In  practical 
units  the  maximum  electron  energy  gain  in  the  ICA  in  terms  of  die  acceleration  gradient  and  pulse 
duration  is 

Un^^JMeV]  <  3  X  10*G(e,  a,  b)E^[MV  /  m]TL[sec] ,  (46) 

where  the  values  of  Ez  and  Xl  are  obtained  from  Eq.  (28).  Using  Eq.  (28)  the  pulse  duration  in 
Eq.  (46)  can  be  written  in  terms  of  the  accelerating  gradient  Hence,  the  maximum  electron 
energy,  given  by  Eq.  (46),  can  be  expressed  in  terms  of  only  the  accelerating  gradient 

Vin.  Results 

We  first  describe  the  results  for  the  ICA  configuration  with  a  dielectric  liner  for  three 
wavelength  drivers,  X  =  10  iim,  1  mm,  and  1  cm.  The  plots  principally  show  the  scaling  of  the 
critical  time  with  the  accelerating  gradient  as  the  parameters  are  varied.  The  minimum  value  of 
tent  shown  in  all  the  figures  corresponds  to  a  pulse  duration  of  30  wave  periods.  Technological 
issues  as  well  as  longitudinal  dispersion  limit  the  pulse  duration  to  several  tens  of  wave  periods. 
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For  a  X  =  10  nm,  1  mm,  and  1  cm  wavelength  driver,  the  minimum  value  of  ‘Cmt  used  in  the  figures 
is  1  psec,  0.1  nsec,  and  1  nsec,  respectively. 

Figure  4  is  a  plot  of  the  critical  lime  as  a  function  of  the  accelerating  gradient,  for  X  =  10 
p.m,  as  the  inner  radius  of  the  waveguide  is  varied.  The  dielectric  layer  thickness  b-a  =  50  pm  in 
all  three  cases,  with  a  =  30, 50,  and  100  pm  for  the  three  curves.  Figure  5  shows  the  variation  of 
as  a  function  of  E*,  for  X  =  10  pm,  for  various  values  of  initial  electron  density  in  the  dielectric 
liner.  The  addition  of  a  small  number  of  electrons  to  seed  the  collisional  ionization  process 
reduces  the  allowable  accelerating  gradient  Er  for  a  given  value  of  tent-  Rgurc  6  shows  the 
variation  of  Tent  versus  the  accelerating  gradient,  for  X  =  10  pm,  for  several  values  of  ionization 
energy  associated  with  the  dielectric  liner  material.  For  a  given  critical  pulse  time  Tmt,  higher 
ionization  potentials  permit  higher  accelerating  gradients.  Figure  7  plots  the  critical  time  as  a 
function  of  the  accelerating  gradient,  for  X  =  10  pm,  as  the  dielectric  constant  e  is  varied.  As  the 
dielectric  constant  is  increased,  the  electric  field  in  die  dielectric  decreases  (to  maintain  continuity 
of  the  displacement  field  £Er  across  the  vacuum-dielectric  boundary),  and  thus  the  ionization  level 
is  reduced.  It  should  be  noted,  however,  that  for  high  values  of  e,  ionization  is  likely  to  occur  at 
die  surface  rather  than  in  the  bulk  of  the  dielectric  layer,  and  thus  the  high  gradients  predicted  for 

large  values  of  e  result  in  surface  breakdown. 

Figure  8  is  a  plot  of  the  maximum  electron  energy  Unua  as  a  function  of  the  accelerating 
gradient  E*  from  Eq.  (46).  In  obtaining  this  plot  the  driver  pulse  duration  Tl  is  set  equal  to  Tmt 
given  by  Fig.  4.  For  energy  gains  above  Umax  ~  10  MeV,  Rg.  8  indicates  that  die  accelerating 
gradient  is  limited  to  E^  -  20  MV/m.  This  rather  low  value  of  accelerating  gradient  corresponds 
to  the  onset  of  collisional  ionization  (see  Fig.  3)  for  these  parameters. 


20 


We  now  consider  the  example  of  a  millimeter  wave  driver.  Rgure  9  shows  the  variation  of 
XaH  vwth  E*.  for  X  s  1  mm,  as  the  radius  of  the  innerwall  radius  is  varied.  The  scaling  here  is 
similar  to  the  X  s  10  |im  case  shown  in  Figure  4,  but  the  allowable  accelerating  gradient  is  lower. 
The  longer  pulse  lengths,  compared  to  the  10  jim  wavelength  example,  are  the  principle  reason 
for  the  lower  accelerating  gradients.  Recall  diat  the  minimum  value  of  Xan,  in  all  the  examples,  is 
taken  to  be  equal  to  30  wave  periods  of  die  driving  field.  This  is  an  arbitrary  but  technologically 
reasonable  lower  limit  on  the  pulse  duration.  Figure  10  is  a  plot  of  w  versus  E»  for  X  =  1  mm,  as 
die  seed  electron  density  is  varied.  This  case  is  similar  to  the  X  =  10  pm  case  shown  in  Fig.  5. 

The  addition  of  a  small  amount  of  seed  electrons  reduces  the  accelerating  gradient  significantly. 

A  Case  1  example  is  shown  in  the  next  figure,  where  the  driver  wavelength  is  larger  than  die  inner 
dielectric  radius,  X  >  a.  Figure  1 1  shows  toit  versus  Ez  for  X  =  1  mm  and  initial  (seed)  electron 
density  npo  =  0  and  1  cm*^  The  lower  limit  on  Ez  is  due  to  the  sudden  turn  on  of  collisional 
ionization  rate  W  as  a  function  of  electric  field,  see  Rg.  3.  Rgure  12  shows  the  maximum 
electron  energy  Unm  as  a  function  of  the  accelerating  gradient  Ez  for  die  millimeter  wave  example 
in  Fig.  9.  Figure  12  indicates  that  for  energy  gains  in  excess  of  Ubk  -  20  MeV  the  accelerating 
gradient  is  less  than  Ez  ~  20  MV/m. 

The  solid  curve  in  Fig.  13  shows  Xait  versus  E*  for  a  driver  wavelength  of  X  =  1  cm  and 
a  =  0.2  cm,  corresponding  to  Case  1.  The  dashed  curve  in  Rg.  13  shows  the  maximum  electron 
energy  Ub«  as  a  function  of  the  accelerating  gradient  E*.  Figure  13  indicates  that  for  energy  gains 
above  Umix  ~  10  MeV  the  accelerating  gradient  is  less  than  Ez  -  100  MV/m,  due  to  the  onset  of 
coUisional  ionization. 
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Figures  14  and  15  show  the  critical  pulse  duration  Xait  versus  the  accelerating  gradient  Ej 
for  a  gas-filled  ICA.’Thegas  is  taken  to  be  helium  having  an  ionization  potentialof  Ui  =  24.6  eV. 
The  variation  of  the  critical  time  with  tiie  accelerating  gradient  Ez,  at  X  -  10  urn,  is  shown  for 

gas  pressures  Pj  of  3, 10,  and  30  atmospheres.  Since  E“l=AEs=7xlO*atl  atm  for  helium,  the 
dispersion  relation  for  the  lowest  order  mode,  po,  =  2.4,  gives  a  =  265, 145,  and  84  jtm 
respectively  for  the  three  pressures.  This  gas-filled  case  has  a  lower  accelerating  gradient  tiian  the 
dielectric  liner  case  shown  in  Rg.  4,  because  for  phase  velocities  equal  to  c,  the  dispersion  relation 
requires  the  wall  radius  to  be  so  large  that  E*  «  E/r  =  a).  In  the  dielectric  liner  cases,  Ez  can  be 
the  larger  field  component  Hgure  15  is  a  plot  of  tail  versus  Ez  for  a  X  =  1  mm  driver  with  the 
same  three  helium  gas  pressures  shown  in  Fig.  14.  The  wall  nulius  a  =  26.5, 14.5,  and  8.4  mm  for 
gas  pressures  Pj  of  3, 10,  and  30  atmospheres,  respectively.  The  allowable  accelerating  gradient 
is  a  few  MV/m.  This  is  due  in  part  to  the  large  a/X  ratio  required  to  have  phase  velocities  equal  to 
c.  In  addition,  the  coUisional  ionization  (avalanche)  time  1/W  even  at  tfiese  modest  fields  is  much 
less  than  the  pulse  duration.  Hgure  16  shows  the  maximum  electron  energy  Ubk  as  a  function  of 
the  accelerating  gradient  Ez  for  the  gas-filled  ICA  using  the  parameters  in  Rg.  14,  wth  gas 
pressures  3, 10,  and  30  atm.  Comparing  Hg.  16  with  Rgs.  8, 12,  and  13  we  find  that  higher 
acceleration  gradients  can  be  obtained  in  the  gas-fiUed  ICA  to  reach  a  given  maximum  energy  as 
compared  to  the  dielectrically-lined  ICA  configuration. 

IX.  Summary  and  Discussion 

The  analysis  in  this  paper  implies  that  the  combination  of  ionization  and  pulse  lethargy 
effects  impose  severe  limitations  on  the  energy  gain  in  inverse  Cherenkov  accelerators.  Two 
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configurations  for  the  ICA  have  been  considered  in  which  the  electromagnetic  driver  (e.g.,  laser, 
millimeter  or  centimeter  wave  source)  is  propagated  in  a  waveguide  that  is,  i)  lined  with  a 
dielectric  material  or  ii)  filled  with  a  neutral  gas.  Using  a  dielectrically-lined  or  gas-filled 
waveguide,  i)  avoids  diffraction  of  the  electromagnetic  driver  beam,  ii)  overcomes  phase  velocity 
slippage,  and  iii)  the  acceleration  is  in  vacuum  (for  the  dielectrically-lined  configuration).  In  either 
configuration  the  electromagnetic  driving  field  has  an  axial  electric  field  with  phase  velocity  equal 
to  c.  The  intensity  of  the  driver  in  the  ICA,  and  therefore  the  acceleration  gradient,  is  limited  by 
tunneling  and  collisional  ionization  effects.  Partial  volume  ionization  of  the  dielectric  liner  or  gas 
can  lead  to  significant  modification  of  the  electromagnetic  properties  of  the  waveguide,  altering 
the  phase  velocity  of  the  accelerating  field  and  causing  phase  velocity  slippage,  thus  disrupting  the 
acceleration  process.  Ionization  effects  therefore  impose  an  upper  limit  on  the  driver  pulse 
duration.  The  driving  pulse,  however,  cannot  be  arbitrarily  short  in  duration  ance  the  pulse, 
having  a  group  velocity  less  than  c,  slips  behind  the  accelerated  electrons.  This  effect  is  referred 
to  as  group  velocity  slippage  or  pulse  lethargy  and  places  a  lower  limit  on  the  pulse  duration. 
Limitations  on  the  driver  pulse  duration,  accelwating  gradient  and  maximum  accelerated  electron 
energy,  due  to  ionization  effects  and  pulse  letfiargy,  have  been  obtained. 

To  obtain  these  limitations  a  dispersion  relation  for  transverse-magnetic  (TM)  modes  of  a 
dielectrically-lined  and  gas-filled  optical  waveguide  has  been  derived.  The  TM  mode  consists  of  a 
radial  and  an  axial  electric  field  as  well  as  an  azimuthal  magnetic  field.  The  anal  electric  field, 
which  is  responsible  for  acceleration,  is  peaked  on  axis,  while  the  radial  electric  field  vanishes  on 
axis.  It  is  shovm  that  when  die  average  plasma  density  reaches  a  critical  value,  due  to  ionization, 
the  optical  properties  of  the  waveguide  are  modified  and  the  phase  velocity  of  the  electromagnetic 
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wave  is  significandy  altered,  lesulting  in  phase  velocity  slippage.  The  buildup  of  the  plasma  within 
the  dielectric  material  is  analyzed  using  a  rate  equation  wduch  includes  tunneling  and  collisional 
ionization.  We  have  obtained  an  exfwession  for  the  limit  on  the  electromagnetic  pulse  duration  in 
terms  of  the  accelerating  gradient  for  both  the  dielectrically-lined  and  gas-filled  waveguide  ICA 
configurations.  To  avoid  the  effects  of  ionization,  the  pulse  duration  of  the  electromagnetic  driver 
%L  in  the  ICA  must  be  less  than  a  critical  pulse  duration  time  Xan  given  by  Eq.  (28).  The  value  of 
XaH  is  a  sensitive  function  of  the  peak  electric  field  and  thus  the  accelerating  gradient  The 
limitation  on  the  pulse  length  due  to  ionization,  Eq.  (28).  together  with  the  limitations  imposed  by 
pulse  lethargy,  Eqs.  (43)  and  (45),  present  a  severe  limit  on  the  interaction  length  and  therefore 
the  mtodmum  energy  gain,  Eq.  (46),  in  the  ICA.  The  maximum  accelerated  electron  energy  can 
be  expressed  in  terms  of  only  the  accelerating  gradient  Our  results  indicate  that  relatively  modest 
accelerating  gradients  can  be  achieved  in  the  ICA. 
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onductor 
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Rg.  1.  Schematic  of  a  dielectrically-lined  optical  waveguide.  The  dielectric  material  with 

dielectric  constant  e  lies  between  the  vacuum  (0  <  r  $  a)  region  and  perfect  conductor 
at  r  =  b. 
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Fig.  2.  Axial  and  radial  components  of  the  electric  field  for  a  dielectrically-lined  waveguide 
supporting  an  axially  symmetric  TM  mode  with  phase  velocity  equal  to  c.  The  solid 
(dotted)  curve  denotes  the  radial  (axial)  component  of  the  electric  field  normalized  to 
the  peak  radial  field  at  r  =  a.  The  parameters  arc  X  =  10  jim,  a  =  30  p.m,  b  -  a  =  50 
^m,  and  £  =  3. 
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Fig.  3.  Plot  of  tunneling  ionization  source  function  S  and  coUisional  ionization  rate  W  as  a 
function  of  normalized  electric  field  E  =  E  /  Ey.  In  this  plot  Ui  =  8  eV,  X  =  10  |ira. 


Vm  =  sec*',  and  Vo  =  lO”  sec*'. 


29 


Ej  (GV/m) 


Fig.  4.  Plot  of  critical  pulse  time  lent  versus  the  accelerating  gradient  Ei  for  a  laser  of 

wavelength  X  =  10  Jira,  inner  radius  a  =  30  pm  (solid  curve),  50  pro  (dotted  curve), 
and  100  pm  (dashed  curve),  b  -  a  =  50  pm.  In  addition,  £  =  3,  Ui  =  8  eV,  npo  =  0,  Vm 
lO'*  sec’'  and  Vo  =  lO"  sec’*. 
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Fig.  6.  Plot  of  critical  pulse  time  Xaii  versus  the  accelerating  gradient  Ej  for  a  laser  of 

A 

wavelength  X  =  10  nm.  The  ionization  potential  Ui  =  4  eV  (dotted  curve),  8  eV  (solid 

curve),  12  eV  (dashed  curve).  The  other  parameters  are  the  same  as  in  Fig.  5  with  npo 

=  0. 
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Fig.  5.  Plot  of  critical  pulse  time  Xcni  versus  the  accelerating  gradient  E*  for  a  laser  of 

wavelength  X  =  10  |Xni.  The  initial  electron  densi^  npo  is  0  (solid  curve),  1  cm 
(dotted  curve),  and  10*  cm'*  (dashed  curve).  In  addition,  a  =  30  pm,  b  =  80  pro,  £  =  3, 
Ui  =  8  eV,  v„  =  lO’*  sec*'  and  Vo  =  lO"  sec*'. 
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Fig.  7. 


Plot  of  critical  pulse  time  Xoh  versus  the  accelerating  gradient  E*  for  a  laser  of 
wavelength  5l  =  10  lira.  Ihe  dielectric  constant  e  =  1.5  (dotted  curve),  3  (solid  curve), 
9  (dashed  curve).  The  other  pararaeters  are  the  same  as  in  Fig.  6  with  Ui  =  8  eV. 
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Ez  (GV/m) 

Fig.  8.  Plot  of  maximum  electron  energy  versus  accelerating  gradient  for  the  parameters  in 
Fig.  4  with  a  =  30  pm  (solid  curve),  50  pm  (dotted  curve),  and  100  pm  (dashed 
curve). 
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Plot  of  critical  pulse  time  Tent  versus  the  accelerating  gradient  Er  for  a  millimeter  wave 
driver,  X  =  1  mm.  The  inner  radius  a  =  3  mm  (solid  curve),  5  mm  (dotted  curve),  and 
10  mm  (dashed  curve),  where  b  -  a  =  5  mm.  In  addition,  £  =  3,  Ui  =  8  eV,  npo  =  0,  Vm  = 
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Plot  of  critical  pulse  time  Tait  versus  the  accelerating  gradient  E*  for  a  millimeter  wave 
driver,  X  =  1  mm.  The  initial  electron  density  iipo  is  0  (solid  curve),  1  cm*^  (dotted 
curve),  and  10*  cm**  (dashed  curve).  The  parameters  are  the  same  as  in  Fig.  9,  with  a 


=  3  mm  and  b  =  8  mm. 


cm  ins 


s-1 


Ej  (GV/m) 

Rg.  1 1.  Plot  of  critical  pulse  time  Xan  versus  the  accelerating  gradient  Ez  for  a  millimeter  wave 

driver,  X  =  1  mm.  The  initial  electron  density  npo  is  0  (solid  curve)  and  1  cm’^  (dashed 
curve).  The  parameters  are  a  =  0.2  mm,  b  =  0.283  mm,  e  =  10,  Ui  =  8  eV,  Vm  =  10*^ 
sec'‘  and  Vo  =  lO"  sec'‘. 
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Plot  of  critical  pulse  time  tent  versus  the  accelerating  gradient  Ex  (solid  curve)  and 
Ubuix  versus  Ex  (dashed  curve)  for  X  =  1  cm  and  a  =  0.2  cm.  In  addition,  b  -  a  = 
0.0833  cm,  e  =  10,  Ui  =  8  eV,  iino  =  0,  v„  =  lO'*  sec*  and  Vo  =  lO’*  sec*. 
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Fig.  14.  Plot  of  criUcal  pulse  time  Tent  versus  the  accelerating  gradient  Ex  in  a  gas  (helium)- 
fiUed  ICA  with  laser  wavelength  X  =  10  nm.  The  gas  pressures  and  corresponding 
wall  radii  are  3  atm  and  265  pm  (solid  curve),  10  atm  and  145  pm  (dotted  curve),  and 
30  atm  and  84  pm  (dashed  curve),  with  Ui  =  24.6  eV,  Vb/P*  =  lO'*  sec'Vatm  and  VofP g 
=  (2m/M)VB/Pg  =  3  X  10*  sec’Vatm. 
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Fig.  15.  Plot  of  critical  pulse  time  'Coii  versus  the  accelerating  gradient  Ez  in  a  gas  (helium)- 

filled  ICA  with  laser  wavelength  X  =  I  mm.  The  gas  pressures  and  corresponding  wall 
radii  are  3  atm  and  26.5  mm  (solid  curve),  10  atm  and  14.5  mm  (dotted  curve),  and  30 
atm  and  8.4  mm  (dashed  curve),  with  Ui  =  24.6  eV,  Vn/P*  =  lO'^  sec'Vatm  and  Vo/Pg  = 
(2m/M)VB/Pj=  3  X  10*  sec'Vatm. 
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Fig.  16.  Plot  of  maximum  electron  energy  versus  accelerating  gradient  for  the  parameter  in  Fig. 

14  with  Pj  =  3  atm  (solid  curve),  10  atm  (dotted  curve),  and  30  atm  (dashed  curve). 
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